Avalanche Dynamics in Wet Granular Materials 
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We have studied the dynamics of avalanching wet granular media in a rotating drum apparatus. 
Quantitative measurements of the flow velocity and the granular flux during avalanches allow us to 
characterize novel avalanche types unique to wet media. We also explore the details of viscoplastic 
flow (observed at the highest liquid contents) in which there are lasting contacts during flow, leading 
to coherence across the entire sample. This coherence leads to a velocity independent flow depth at 
high rotation rates and novel robust pattern formation in the granular surface. 

PACS numbers: 45.70.-n, 45.70.Ht, 45.70.Mg 



Avalanches and landslides [|l|, |2|, |^ are among the most 
dramatic of natural catastrophes, and they also provide 
an evocative metaphor for a wide range of propagating 
breakdown phenomena from impact ionization in semi- 
conductors to magnetic vortex motion in superconduc- 
tors [||. On the other hand, the existence of avalanches, 
i.e. the sudden collapse of the system previously frozen 
into a high energy state, is a fundamental manifestation 
of the metastable nature of granular materials. Studies of 
avalanches [^, 0, || and surface flows |9| in granular media 
have largely focused on dry grains. By wetting such me- 
dia, however, one introduces controllable adhesive forces 
between the grain s |iq|, 11| w hich lead to qualitatively 
new behavior |13|, y, |15|, |6[ In our previ- 
ous work we identified three fundamental regimes for 
the repose angle of wet granular materials as a function 
of the liquid content. The granular regime at very low 
liquid contents is dominated by the motion of individual 
grains; in the correlated regime corresponding to inter- 
mediate liquid contents, a rough surface is formed by the 
flow of separated clumps; and the repose angle of very 
wet samples results from cohesive flow with viscoplastic 
properties. 

Here we report investigations of the avalanche dynam- 
ics and flow properties of wet granular materials, em- 
ploying a rotating drum apparatus (a cylindrical cham- 
ber partly filled with a granular medium and rotated 
around a horizontal axis) [^,0. At low rotation rates, 
the medium remains at rest relative to the drum while 
its surface angle is slowly increased by rotation, up to a 
critical angle {Omax) where an avalanche occurs, thus de- 
creasing the surface angle to the repose angle {6r). The 
flow becomes continuous at high rotation rates, but the 
transition between avalanching and continuous flow is 
hysteretic in rotation rate in dry media [||, ^ . 

Previous studies of cohesive granular media in a rotat- 
ing drum [l^, Q have focused on the surface angles 
of the medium before and after avalanches. For example, 
Quintanilla et al. recently performed a statistical analy- 
sis of avalanche size based on these angles, and showed 



that the average clump size increased with cohesion [gO| . 
In our measurements, we focus instead on characteriz- 
ing the dynamics of cohesive flow. We quantitatively 
investigate the flow dynamics during avalanches at dif- 
ferent liquid contents by analyzing the time evolution 
of the averaged surface profile obtained from hundreds 
of avalanche events, and we also measure surface veloc- 
ities during continuous flow. In particular, we explore 
the nature of the viscoplastic flow, which displays unique 
characteristics associated with coherent motion over the 
entire granular surface. 

We studied glass spheres thoroughly mixed with small 
quantities of hydrocarbon oil |2^] varying between r = 
0.001% and 5% of the void volume. In this regime the 
flow of oil due to gravity can be neglected. The inner 
diameter of our drum was 16.8 cm, its width was 3.2 
cm, and the granular filling was 30%. Measurements 
performed in a thinner (2 cm) drum verified that wall 
effects do not modify the qualitative behavior. We stud- 
ied two different sizes of beads {di = 0.5mm±20% and 
d2 = 0.9mm±ll%), and the results were qualitatively 
equivalent in the two cases except where noted otherwise. 
The experiment was recorded by a CCD video camera 
interfaced to a computer which could analyze the spa- 
tiotemporal evolution of the surface profile (height varia- 
tions in the axial direction were negligible). To determine 
the position of the surface, we used background illumina- 
tion and an algorithm based on local contrast which gave 
us a resolution of 5t = 0.03 seconds and 5x = 0.5 mm. 

Surface angles. — We first investigated the avalanching 
- continuous flow transition for various liquid contents. 
We slowly increased and then decreased the rotation rate 
and measured the average surface angles, Omax and 6r in 
the avalanching regime. We consider the flow continuous 
when the medium never comes to rest with respect to the 
drum, which coincides with = Omax ~ 

l9. < 1°. The 

results of typical runs are presented in Fig. [| (a), (b) 
and (c), where Omax and Or are plotted as a function of 
the rotation rate. 

The three basic types of behavior observed previously 
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FIG. 1: (a) (b) (c) Omax and Or as a function of rotation rate 
for 3 different liquid contents. The error bars represent the 
standard deviation of the observed angle distributions, (d) 
The maximum angle, Omax, the average avalanche size, (A^), 
and the width of the avalanche size distribution 2aAe as a 
function of liquid content. (The data were taken using small 
beads, d — 0.5 mm) 



are also reflected in the present experiment. In the gran- 
ular regime at low liquid contents, (e.g. r = 0.04%), the 
behavior is qualitatively similar to the dry case, and we 
observe a clear hysteretic transition between continuous 
flow and avalanching (Fig. |l| (a)). At a somewhat higher 
liquid content (r = 0.125%, see Fig. ^ (b)) the cor- 
related behavior is marked by a lack of hysteresis, and 
avalanching behavior is observed at a relatively high ro- 
tation rate compared to the other regimes. The continu- 
ous flow in this regime consists of a stream of separated 
clumps rather than the constant flux seen in the other 
two regimes. At r = 0.375% (Fig. |l| (c)) hysteresis is 
again observed, reflecting the onset of the viscoplastic 
continuous flow, which is smooth and coherent over the 
entire sample. 

An interesting feature of the curves in Fig. (a-c), 
is that Omax decreases with increasing rotation rate, al- 
though at slow rotation (where the rotation of the drum 
during the avalanche is negligible) the rotation rate influ- 
ences the avalanche process only through changes in the 
waiting time between successive avalanches. We find that 
Omax depends logarithmically on this waiting time [25|, 
reminiscent of other granular aging phenomena 22 1. 
The effect is more pronounced at higher liquid contents, 
and we attribute it to the motion of liquid flowing to- 
wards the contact points rather than condensation effects 
[p!5| since the vapor pressure of the oil is low. 



If the rotation rate is sufficiently low, then discrete 
avalanches occur at all liquid contents. In the fol- 
lowing paragraphs we analyze this avalanching regime. 
Fig. (d) shows {Omax) (averaged over several hundred 
avalanches) as a function of liquid content. The distribu- 
tion is close to Gaussian, and the error bars in the figure 
indicate the standard deviation. The average avalanche 
size (A^) = {Omax — Or) and the width of the avalanche 
size distribution 2crA^ = 2((A6>2) - (Al9)2)^/2 are also 
plotted in Fig. |l| (d). We interpret these results be- 
low in the context of a more detailed investigation of the 
avalanche dynamics. 

Flow Dynamics. — By using the rotating drum ap- 
paratus, we can obtain information not just about the 
medium before and after the avalanche events, but we 
can also study the details of the grain motion dur- 
ing avalanche events. In order to analyze the dynam- 
ics of avalanches we have obtained two dimensional 
space-time matrices, h{x^t)^ characterizing the sample 
surfaces throughout the avalanche process which can 
then be analyzed to produce a variety of information 
about the individual avalanches. By taking deriva- 
tives of the h{x^t) data, we obtain the local angle 
{a{x,t) = arctan[9a^/i(x, t)]) and the local vertical ve- 
locity {u{x,t) = dth{x,t)) of the surface proffie. Fur- 
thermore, by then integrating the vertical velocity (using 
the continuity equation and assuming constant density 
[p3[ ), we also obtain the local flux in the avalanche, i.e. 
(j){x^ t) = pw J^jj^2 9th{x' ^ t)dx' (where p is the grain den- 
sity, w is the width of the drum, and D is the drum di- 
ameter), which represents the material flowing through a 
given vertical plane at position x. 

In Fig. ^ (a) we present snapshots of the progression 
of single avalanches for several typical liquid contents. 
Fig. H (b) displays a{x^t) as a function of space (hori- 
zontally) and time (downwards) for the same individual 
avalanches. In Fig. || (c-e), we present the average be- 
havior of 300-500 avalanches at the same liquid contents, 
and show similar graphs of the time evolution of (a(x, t)), 
(ii(x,t)), and ((/)(x,t)), where () denotes averaging over 
avalanches. By obtaining these quantitative measures 
of the averaged properties, we can separate the robust 
characteristics of the avalanche dynamics from the large 
fluctuations which are inherent in avalanche processes. 

In the avalanches among dry grains [^, |^, ||], the 
avalanches consist of flow in a thin layer near the smooth 
surface and have a much longer duration and much 
smaller flux than in the wet media. The size of the 
avalanches is A^ 2° ± 0.5°. At low but nonzero liquid 
contents (e.g. r = 0.04%) Omax and A^ become larger 
due to the onset of intergrain cohesion. In this granu- 
lar regime the avalanche is initiated by a front traveling 
downhill (bright line in Fig. || (c)), but we can also ob- 
serve a second front traveling uphill (dark line). The 
second front is also apparent in the motion of the point 
of maximum flux. Fig. (e), and it corresponds to the 



3 



Dry 



T = 0.04% T = 0.12% T = 5.00% 



m 



J § 




FIG. 2: Dynamics of avalanches of different types with grain 
size d — 0.5 mm. (a) Snapshots (at 0.1—0.2 s time intervals) of 
four single avalanches corresponding to different liquid con- 
tents. The third dimension is used for a brightness-coded 
representation of the local slope (a). The green lines indicate 
approximate slip planes in the correlated regime, and the red 
line shows the traveling quasi-periodic surface features in the 
viscoplastic regime, (b) The local slope with the same bright- 
ness coding, as a function of space and time. A horizontal line 
corresponds to a surface profile at a given instant, and time 
propagates downwards (thus avalanches propagate leftwards 
and down). The stripes indicate lasting surface features, (c) 
(d) (e) The characteristic features of the avalanches averaged 
over several hundred independent avalanches. The displayed 
quantities are: (c) local surface angle, (d) the rate of change 
of local height, and (e) the local grain flux (the white lines 
indicate the point of maximum current). Note that the sur- 
face patterns at the highest Uquid contents are robust against 
averaging. The standard deviations of the averaged data sets 
were typically 5 — 25 % of the mean value. In some cases, 
however, the noise originating from taking numerical deriva- 
tives leads to standard deviations comparable to the mean. 
For this reason averaging over several hundred independent 
events was absolutely necessary, resulting in standard errors 
ranging from 0.2 — 6 %. 



grains reaching a solid barrier at the bottom [ p4[ . 

At intermediate liquid contents (the correlated regime, 
0.1% ^ T < 0.5%), both dmax and increase sub- 
stantially. In this regime, the principal failure mecha- 
nism is a fracture along a curved slip plane (approxi- 
mated by the green lines in Fig. || (a)), analogous to 
the dynamics of a class of geological events known as 
"slides" [||. In the figure there is a single slip plane, 
which leads to large avalanches with a relatively nar- 
row distribution. However, at larger liquid contents the 
avalanches occur through a succession of local slip events 
followed by a large avalanche extending over the whole 
surface This leads to large fluctuations in the size 

of the avalanches (see Fig. ^ (d)) as observed in cohe- 
sive powders The medium becomes more cohesive 
with increasing liquid content, and the upward traveling 
front disappears for r > 0.2% since the material moving 
as a connected block stops coherently when it hits the 
bottom. 

At the highest liquid contents, r 0.5 — 1%, the onset 
of coherent viscoplastic flow is apparent through numer- 
ous qualitative changes in the dynamics of flow. The co- 
herence of avalanche flow in this regime is demonstrated 
by the parallel lines during the avalanche in Fig || (b) 
and (c). This behavior is qualitatively similar to a differ- 
ent class of geological event, called a "debris flow" or 
"mudflow" H, |[. Since the whole surface moves co- 
herently (rather than breaking apart and evolving sep- 
arately in different parts of the surface), fluctuations are 
strongly suppressed Q and the avalanche size distribu- 
tion is rather narrow (Fig. ^ (d)). The small decrease in 
^max is probably due to lubrication jl^. 

A novel property of the viscoplastic avalanches is the 
topology of the top surface which spontaneously forms 
a nearly periodic pattern (seen in Fig. || (a) and (b)). 
This surface structure is maintained essentially intact 
during the avalanche, indicating that there are lasting 
contacts in the flowing layer. Moreover, the same pat- 
tern in the surface features is reproduced at the end of 
each avalanche. This is demonstrated most clearly in Fig 
^ (c), where the average of 347 avalanches of the r = 5% 
sample has the same features as the typical individual 
avalanche shown in Fie;. I (a) and (b). The robust na- 
ture of the surface structure of the wettest grains is in 
sharp contrast to the other regimes where averaging com- 
pletely smoothes out the smaller surface features. We 
can understand this behavior as resulting from coher- 
ence of the entire flow, which strongly reduces fluctu- 
ations in this regime. With minimal fluctuations, the 
final surface structure after each avalanche is essentially 
the same, thus setting the same initial condition for the 
next avalanche. With the same initial conditions for each 
avalanche, naturally the surface features are reproduced 
each time. Our experiments with the larger (d = 0.9 mm) 
beads also revealed pattern formation, but with a smaller 
characteristic size corresponding to 8-10 grain diameters. 
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FIG. 3: The surface velocity v during continuous flow as 
a function of rotation rate u for dry grains and in the vis- 
coplastic regime, r = 5% (large beads, d = 0.9mm). The 
continuous lines represent power-low fits to the data. The 
linear behavior for the viscoplastic flow indicates that flow 
depth is independent of the rotation rate. 



The difference is probably due to the smaller ratio of the 
cohesive forces to the gravitational forces on the grains. 

In addition to differences in the avalanching behav- 
ior, the viscoplastic regime also displays continuous flow 
which is rather different from that at lower liquid con- 
tents. We compare the viscoplastic and granular con- 
tinuous flows by measuring their surface velocities v (by 
means of tracer particles) as a function of rotation rate 
in Fig. ^ (continuous flow is only observed in the corre- 
lated regime above 6 rpm, and we therefore do not include 
data from that regime here). For these measurements we 
used the larger, {d = 0.9 mm) beads. The viscoplastic 
flow is slower than the granular flow by a factor of two, 
and the two curves correspond to rather different func- 
tional forms. The curvature in the granular case suggests 
that the flow depth increases with rotation rate as pre- 
viously suggested [|2^. The linear v{uj) function of the 
viscoplastic flow suggests a constant flow depth which is 
independent of the flow rate. This observation is consis- 
tent with the coherent nature of the flow: the flow depth 
is not determined by local mechanisms, but is fixed by 
the geometry of the whole system Q . 

Discussion. — While the flows we observe appear 
to have analogies with geological events, it is impor- 
tant to note that real geological materials usually con- 
sist of polydisperse irregular particles often with very 
high (r ^ 100%) water contents and that the scaling 
of our system to geological lengthscales is non-trivial. 
Furthermore, avalanche studies in real soil have demon- 
strated additional phenomena associated with soil sat- 
uration [p^. Interestingly we still recover some of the 



basic dynamical processes in our model system, which 
constitutes an important step towards the description of 
qualitatively different flow behaviors in the framework 
of a single model |2^. The changes in the dynam- 
ical behavior with wetting are associated with the in- 
creasingly coherent nature of the flow, i.e. the forma- 
tion of coherently moving clusters - clumps - due to the 
increased cohesion and viscous effects. Within a clus- 
ter, local velocity fluctuations should be suppressed, and 
thus the local granular temperature (T = {v'^) — (v)'^) 
should approach zero, but the clusters themselves both 
form and break apart during an avalanche process in a fi- 
nite container. An important theoretical question raised 
by our data is how a length scale describing the size of 
the clumps may emerge from a granular flow model, and 
how such a length scale would vary with the type of me- 
dia, the nature of intergranular adhesion, and the type 
of granular flow. 
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